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Abstract 

The calculation of heat-kernel coefficients with the classical DeWitt algo- 
rithm has been discussed. We present the explicit form of the coefficients up 
to h$ in the general case and up to /i™ n for the minimal parts. The results 
are compared with the expressions in other papers. A method to optimize 
the usage of memory for working with large expressions on universal computer 
algebra systems has been proposed. 
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Introduction 



The calculation of the quark determinant modulus can be reduced to the calculation 



where K(r) = e~ Ar is the so-called "heat kernel" for operator A = D^D, with D = 
7^0^ being Dirac differential operator in the presence of the external background 
fields. The coefficients hk of the expansion of the interaction part of the heat kernel 
in powers of the proper time r are known as the heat-kernel coefficients. They 
determine the low-energy expansion of the effective one-loop action [[|. Therefore, 
the calculation of these coefficients is an alternative method to the direct calculation 
of Feynman diagrams 0], thus having a fundamental character. 

This paper is an extension of our previous paper We discuss the implemen- 
tation of the classical DeWitt algorithm || to calculate the higher-order heat-kernel 
coefficients by means of computer algebra. After presenting the results obtained 
with this method we briefly compare them with those of other groups using different 
techniques. This method is demonstrated to obtain also the nonlocal corrections 
taking into account the finite sizes of mesons, that leads to the modification of the 
heat-kernel equation. We also discuss technical questions connected with the work 
with large expressions in Lisp-based computer algebra systems such as Reduce. Our 
techniques allow us to push the calculations of the heat-kernel coefficients to higher 
orders. 

1 General outline of the heat-kernel techniques 

The logarithm of the determinant of a positive definite operator A is defined in 
proper-time regularization with the integral relation Eq. ([!]). The trace Tr is to be 
understood as a space-time integration and a "normal" trace over Dirac, color and 



of the Schwinger proper-time integral [|TJ: 




(1) 



flavor indices. The operator A is an elliptic one in usual cases: 

A = D f D = d^ + a(x) + fi 2 , 

where = + T M , operator describes vector gauge fields, a(x) is a local oper- 
ator without free derivatives. The explicit form of these operators depends on the 
particular model. The effective parameter /z fixes the regularization in the region of 
low momenta. 

The heat kernel K(r) = e~ Ar satisfies the equation 

^K(r) + AK(r) = 

with the boundary condition 

K(t = 0) = 1. 

The asymptotic behavior of A at short distances is defined by the "free" part 
Using the ansatz 

K = K H 

it is convenient to separate from the heat-kernel its "free" part K , 

1 f (x-y) 2 ] 

< x\K \y >=< x\ exp(-(9^ + fi 2 )r)\y >= — - exp |-// 2 r + J' \ , 

which satisfies equation 

^ + A o^o = 0, 

with the boundary condition 

K (r = 0) = 1. 

The "interaction" part H of the heat kernel satisfies the equation 
d 1 

(— + -z^d? + dPdp + a)H(x, y; r) = 0, 
or r 

H(x,y = x;r = 0) = l, (2) 

with z^ = x^ — y^. The differential operator d^ acts on x only. Using now an 
expansion for H{r) in powers of r 

oo 

H(x,y;r) = ^h k (x,y) ■ r k 

k=0 



a recursive relation is obtained from Eq. (0): 



(n + Zpd^hnfay) 



(3) 




log | det (D\ = ~ log(det D^D) = ~\j^y 2 E - 




where T(a,x) = f£° dt e t t a 1 is the incomplete gamma function and = hk(x, y = x). 

2 Recursive algorithm 

In order to construct a recursive algorithm let us use the basic property of the 
heat-kernel coefficients providing the recursive relation (Q) between h n (x,y) and 
h n _i(x,y). The boundary condition for h (x,y) is determined from Eq. (|3]) with 
n = and h-i(x,y) = 0. The task is to find the coincidence limit h n = h n (x,y)\ z=0 . 
We cannot set z = in Eq. (|3|) directly since when acting by differential oper- 
ator d a on this equation we get a nonvanishing contribution: d a (z^d^h n )\ z=0 = 
Qa^i ' d^h n \ z _Q d a h n \ z _Q. 

One can see that the usage of Eq. (H) for the calculation of h n will produce terms 
like d a d/3 . . . h n (x, y)\ z=0 - In order to get the recursive relation for them we apply the 
product of m differential operators d a dp . . . d w on Eq. (|3|) and take limit z = 0: 




dad/s ■■■d u > h n (x,y)\ 



z=0 



in 



1 




n + m 



where (n + m) > and 




d a df3 ...d u -z^\ z=Q - mdad/3 ...d t 



in 



For P. 



we have the recursive relation 




(6) 



with the boundary condition P = 0, where Rp.,.^ a = [dp . . . d^, d a ]. Performing a 
commutation of the differential operator d a successively through the others dp, 
d w , one can move it to the left side. Then the two products with m differential 
operators cancel each other, and only terms with (m — 2) differential operators are 
left over. Finally we will get the recursive relation 

Rpy...<j}; a, ^ /3a ' d^ . . . d^ -\- dp • Ry...u>;a (7) 

with the boundary condition R. a = 0. Thus, one has to use the recursive relation 
(|) starting from m = to calculate h n (x,y)\ z=Q . After each iteration it is necessary 
to commute all differential operators arising from g^gP or P a p... u to the right side up 
to hk(x,y) introducing commutators of the type 

Sfj, [d^, o], S^i U [d/jii $v\j ^a/jiv \da,i •S'/Lti']) etc. 

r — [dfj,, 0?;/], K a ^ u = [rf Q ,r^ tl/ ], Kp a ^ u = [dp, Kain/\ i etc. (8) 

n and m change under these iterations in the following way: either (n — > n — 1 ) , or 
(m — > m — 2), or (n — > n — 1; m — > m + 2). It is easy to show that after 2n iterations 
only h (x,y) remains without differential operators. At the end one gets the desired 
result by setting the limit z = where ho(x,y)\ z=0 = 1. 

Following the strategy outlined above, the calculation of the heat-kernel coeffi- 
cients is straightforward but cumbersome. The lengthy calculations can be performed 
only with computer support. The calculation of the heat-coefficients is a recursive 
process which can be done by computer algebra very conveniently using the recursive 
relations (^-0) until substitutions cannot be made any further. 



3 Heat-kernel coefficients 

The expressions for the heat-kernel coefficients contain a large number of terms which 
can be related to each other with the following transformations: cyclic properties of 
trace, commutator identities removing physically redundant total derivatives, 
renaming of the dummy indices and Jacobi identities for symbols K aVLV . 

It is necessary to reduce the expressions to some minimal basis of linearly inde- 
pendent terms for the final presentation of the results and their comparison with the 
results of other papers. As soon as we know, it is still not possible for nontrivial cases 
to bring all expressions to a canonical form using some identities. The solution of this 
problem is very important since sometimes it is easier to perform a calculation in a 



noncanonical basis than to bring it to a canonical form or test its equivalence with the 
similar expression in another presentation. For example, we had to perform by hand 
the reduction of the effective chiral lagrangians of the order p 6 in the momentum 
expansion to a minimal basis 0, after most of the initial "trivial" simplifications 
had been done automatically with computer algebra. 

So far we do not know an efficient general algorithm which allows to transform the 
total heat-kernel coefficients into some unique minimal basis. For the reduction of the 
heat-kernel coefficients to the minimal basis let us use what is always possible to do 
- test whether two expressions are equivalent or not. Here it is sufficient to express 
all Sfj,.., and V ^ according to their definitions in terms of operators a and d^. The 
comparison of such expressions is easier since the only equivalent transformations 
here are cyclic permutations and renaming of dummy indices. If we can test the 
equivalence of any expressions with some unknown coefficients, then we can expand 
any expression in the minimal basis solving the system of linear equations for the 
coefficients of this expansion. The minimal basis itself can also be constructed by 
rejecting the linearly dependent terms and including the linearly independent ones. 
These operations were fully automated by means of the computer algebra and they 
have allowed us to obtain the heat-kernel coefficients hi, ... h^ in the minimal-basis 
form: 

ho = 1, 
hi = —a, 
Trh 2 = ±Tr | a 2 + ^} , 

Tr h-\ = - Tr ( — a 3 H — Sf, aTf,,. H K?„.„ Y ua T aL \ , 

° g I 2 2 10 15 J 



Tr h A = 1 Tr {a 4 + a 2 S,, + (|a 2 I%, + 1 (aT, 

+ ^ U Y vp^ pp + (^~~^ a -Kuup ~l~ Y§^ aa ^ v ~ YE ^ 1°^ aP^J 
/ 17 2 1 

_|_ f F 2 F 2 _| F F F F _| F F F F 
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Tr k5 = 2i0 Tr { " ^ + (" 2a3 ^ - a2S D + {- a " T l» ~ \ a2T P" aT L ~ 



/2 2 2 \ / 2 2 

+ ( 77° S^K VU ^ — —S^S v aY^ u 1 + I —&S n^Svv — —S^S V 



g •« " p^^-vvp i^p^v^^pv I I I ^^pp^vv g' 

2 
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K^ v K aav 4- ^dK ^^ u QjK aau ^ciS ua 4- i 
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4 4 2 2 2 

~^^S ua Y a pY p v 2T^ ^ jU ^~ l 10 ^ fiv^- otfi ^ fi-^uupJ- ^ a/3 vvo^- /3a 

2 \ / 1 2 3 2 

^^^r^i^o^r^ J 4- I —QiK ^ vua K ^pp a + — Q,K ixvi/aKa.pPn 2g ^""^"z^ 

^2 *^ mi-^vi/t&K- pj3a ' ^,^K uaoi pY u p 4- — S ^K uuaa pT^ 
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37 1 1 
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189 " p«P M / g^ /•* / IH 

^ 1~89 \ 42 ^ ~ Hnvua-Kafi'y^- ^ ft"/ gg K ^ vva K ^jjp^Y a ry 

gg K ^ vva K a pj3~jY ^ gg ^ jxvvo. -^/3/3/i -^-770: ^ "I - -^2Q ^ [ijivvcx^- PPyya J' "I - h.C. 

Here Hermitian conjugation h.c. is defined by 



t = a, (V-) f = 



a 



rr i f = -r fir ^ = -/r 



jti...> 



CtflV.. 



Our expressions for the heat-kernel coefficients ^4 5 proved to be equivalent with the 
results of Refs. |7|, ||. On the other hand, we have found deviations from Ref. in 



the part containing terms with 6 and 8 indices: 



(24Tr/ i4 ) t hisRef.-(24Tr/ i4 ) Re f. 



cid^d I/dwells -\- ad i^df/djfxdi/dfx cid^d^d^d^ ~\~ ad^di/d^d^d^ 

ID ID .LD ID 

16 2 16 

-\- -^-^cid^di/d^di/d^ -^-^adudi/dotdi/dotd^ 
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-^^d 2 ^d u d a d u dpd a dp + -^^d 2 ^d u d a d 2 pd u d a + -^^d^dydadpdydadp 
16 .o 16 .o 16 



^d^dvdadpdvdpda — —d^dvdadpdadvdp + -^-^d^dvdadp-dadpdv 
32 64 32 



:d^d u d^d u d a dpd a dp + -—zd^d u d^d a d u dpd a dp — —7zd^d u d^d a d^d u d a dij 
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32 

+ ^^d^d u d^d a dpd ll dpd a 

These deviations arise from the terms of the expression from Ref. [0] 

2 2 4 

where {A | B \ C} = ABC + CBA. The corresponding terms of our expression are 

(+^s' a( ,r (JJ , + Yo^-^/zaa^r^rp^) . 

We have performed the following checks of our results. First of all, we have used 
the identities <9( Tr h n ) / da = —h n -\ @. Secondly, substituting the explicit form of 
operators a and V ^ in the Nambu-Jona-Lasinio (NJL) model || and evaluating trace 
over Dirac gamma-matrices we derived the effective chiral p 4 - and p 6 -lagrangians 
describing the low-energy processes with mesons and photons ||. In particular, we 
reproduce the structure coefficients Lj of the p 4 -lagrangians obtained in Ref. [[TTJ. We 
also reproduce the effective Euler-Heisenberg lagrangian |TT[ describing the photon- 



photon scattering for particles of both spin and spin 1/2 1 12 



In the case V ^ = we can also present the next orders of the "minimal" parts of 
the heat-kernel coefficients: 

Tr hf n = ^ Tr j - a 5 + 3a 2 Sj + 2aS,aS fl - aS% - b -S,S v S, u + 

Tr K m = ^ Tr ja 6 + (-4a 3 ^ - 6a 2 S>Sj + (y ^ + 

26 18 26 9 17 \ 



Tr ht n = — V- Tr I - a 7 + ( 5a 4 S 2 + 8a 3 S>£„ + -a 2 SyS, 
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1 



The expressions for the heat-kernel coefficients up to Tr /i™ m have been also 
presented in Ref. [fL3|J , and up to Tr/t™ n in Ref. [13]. The minimal coefficients h™ m 
can be easily transformed into a unique minimal basis [|H| : it is necessary to move all 
identical indices in each symbol, e.g. S^^..., to the rest symbols in every term because 
any addition of some total derivative to the lagrangian does not alter its physical 
content. This way we have checked our expressions for the minimal coefficients to 
coincide with other papers. 



4 Nonlocal corrections for the effective lagrangian 



The usage of the recursive relations allows us to extend the local NJL model || by 
taking into account the finite sizes of mesons [15 in a simple way. In this case the 
modulus squared of Dirac operator D receives some additional contributions propor- 
tional to a small parameter a/A 2 characterizing the size of the nonlocal corrections: 



D+D 



(3d 2 + /i 2 + 2T^ + Tj + a 



+ 



a 
A 2 



b + Q a d a + (a + c)d 2 + 2(r M <9 2 + <9 a rvcP)^ 



2 ■ 



where (3 = 1 + 2a/i 2 /A 2 and a, b, c and are local operators (their explicit form 
can be found in Ref. [|15fl ). 



The modified recursive equation for the heat-kernel coefficients h n (for simplicity 
we study only the case without vector fields) reads 



a 
4A 1 



z 2 ch n+1 (x,y) 
\n + z^ + ^ 



2(a + c)(l + X -z^) - 2(i 2 z^ + X -z^a + 



h n {x,y) 



+ a + 8 2 + 



a 



b + (fy + 2Q")0" + (a + c)& 



h n -i(x,y) = . 



Analogously to the algorithm described above, one can obtain the recursive relations 
for d a dp . . . h n (x, y)\ z=0 , where each iteration gives either (n — > n — 1), or (m — > 
m — 1), or (m — > m — 2), or (n — > n — 1; m — > m + 1), or {n — > n — 1; m — > m + 2). It 
is easy to see that functional F = 3n + m reduces at least by 1 after each iteration. 
Therefore, the final result is obtained after 3n iterations [O . 



5 Implementation of calculation in computer al- 
gebra system Reduce 



We have used the Computer Algebra System (CAS) Reduce |Tj| extended by our 
package for the calculations in chiral meson theories |l7| . Reduce suits well for our 
purposes since it is a universal and open system which can be easily extended by the 
users. It is widely used for our calculations of amplitudes of various meson processes 
and for the derivation and transformation of lagrangians from bosonization of the 
NJL model. It allows us to complete the general mathematical and algorithmic base 
of Reduce by the specific data types and operations for this field. In general it has 
provided a convenient computational environment for studying chiral meson theories. 
A more detailed description of our package can be found in Ref. |17] . 



For the calculations described in this paper the following package features are the 
most essential: 

• transformation by cyclic permutation of noncommutative operator products 
under the trace operation, e.g. Tr [ABC) = Tr (BCA); 

• transformation by redefinition of dummy indices, e.g. S^S^Sv = S vtl S v S^ 

• BTgX output of large expressions (on the basis of RLFI package from Reduce 



library 



For the work with the heat-kernel coefficients the following operations with indices 
have been additionally introduced: 



• ordering indices in a symbol: 




• moving indices from the symbols with too many of them: 





• moving apart identical indices to different symbols: 



Tr (S^S UU ) — > Tr (S^S^) + commutator terms. 



As it was already noted, these procedures are not sufficient to transform expressions 
to some unique form. Nevertheless, they perform the overwhelming part of the work 
to reduce the original expressions. 

Universal Lisp-based CAS such as Reduce cannot work with large expressions 
which do not fit into available computer memory. At the same time the intermediate 
expressions can grow very much during the large-scale recursive calculations. This 
problem is effectively tackled by the CAS Form [pl| which works with the terms 
of large expressions "locally" keeping only a relatively small part of the terms in 
the memory and storing the others on the disk. Another alternative is writing a 
specialized program to calculate the heat-kernel coefficients in the language C |2(J . 

We have implemented the following simple method to resolve the problem of 
uncontrolled growth of expressions causing the memory overflow in Reduce: During 
the work of algorithm the expression is regularly (after each iteration) looked through 
and a limit is imposed on the number of terms to be substituted at the next iteration. 
This can be done by a simple change of all differential operators d^, except the first 
n max ones, to some new operator for which no substitutions are set. The parameter 
n max is chosen to provide an optimal loading of memory without its overflow. This 
method has turned out to be effective. It allowed us to calculate the next orders of 
the heat-kernel coefficients. 

Conclusion 

The usage of the classical DeWitt algorithm allows to generate effective lagrangians 
for a wide range of problems in a simple way. In our previous papers || we obtained 
an effective chiral lagrangian from bosonization of the NJL model at 0(p e ) order 



using this method. It also allows to study the nonlocal corrections for the p 4 chiral 



lagrangian fll5| . The usage of the CAS Reduce extended by a specialized package al- 
lowed us to obtain the higher-order heat-kernel coefficients, reduce these expressions 
to minimal basis and compare the results with some other papers. 
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